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A novel type of one-dimensional (ID) photonic crystal formed by the array of periodically located 
stacks of alternating graphene and dielectric stripes embedded into a background dielectric medium is 
proposed. The wave equation for the electromagnetic wave propagating in such a structure is solved 
in the framework of the Kronig-Penney model. The frequency band structure of ID graphene-based 
photonic crystal is obtained analytically as a function of the filling factor and the thickness of the 
dielectric between graphene stripes. The photonic frequency corresponding to the electromagnetic 
wave localized by the defect of photonic crystal formed by the extra dielectric placed on the place 
of the stack of alternating graphene and dielectric stripes is obtained. 

PACS numbers: 42.70.Qs, 78.67.Wj, 78.67.-n, 78.67.Pt 

I. INTRODUCTION 

Photonic crystals are formed by structures with the dielectric constant periodically varying in space Electro- 
magnetic waves in photonic crystals have a band spectrum and a coordinate dependence caused by this periodicity 
of the dielectric constant. Interestingly enough that the solution of Maxwell's equations with a periodic dielectric 
constant, resulting in photonic band-gap structures, is similar to the solution to Schrodinger's equation for a periodic 
potential, resulting in the electron energy band-gap structures in solids. Electromagnetic waves penetrate in the pho- 
tonic crystal similar to the Bloch waves of electrons in a regular crystal. The width of the photonic band gap depends 
on the geometrical parameters of the photonic crystal and the contrast of the dielectric constants of the constituent 
elements [2, 3]. Different materials have been used for the corresponding constituent elements including dielectrics, 
semiconductors and metals @, H, . The photonic crystals with superconducting elements have been studied in 
Refs. [IM3 .The one-dimensional photonic crystals formed by semiconductors were analyzed in Ref. [17]. It is well 
known that the properties of photonic crystals provide an opportunity to manipulate the emission, propagation and 
distribution of light and photonic crystals can be used as frequency filters. The properties of photonic crystals 

were reviewed in Ref. |18| . 

Photonic crystals are different from the regular solid crystals in the following way. While Schrodinger's equation 
describes the regular solid crystals via the scalar wave function, Maxwell's equations for photonic band-gap crystals 
describes the electric or magnetic field, which is a vector corresponding to the transverse electromagnetic waves. 

A novel type of 2D electron system was experimentally obtained in grap hene, which is a 2D honeycomb lattice of 
the carbon atoms that form the basic planar structure in graphite [19l . |20| . Due to unusual properties of the band 
structure, electronic properties of graphene became the object of many recent experimental and theoretical studies 
[l9l - l25j . Graphene is a gapless semiconductor with massless electrons and holes which have been described as Dirac- 
fermions [26j . The unique electronic properties of graphene in a magnetic field have been studied recently (27l - [30| . 
The space-time dispersion of graphene conductivity was analyzed in Ref. [3l|. A graphene-based two-dimensional 
photonic crystal was proposed in Ref. 1321 . where its frequency band structure was studied. Different types of 2D 
photonic crystals were reviewed in Ref. [331 ] . 

In this Paper, we consider a one-dimensional photonic crystal formed by an array of periodically located parallel 
stacks of alternating graphene and dielectric stripes embedding into a background dielectric medium. The graphene 
stripes are placed one under the other with the dielectric stripes placed between them. We calculate the frequency 
band structure of such a photonic crystal. The photonic band structure usually can be obtained using numerical 
calculations. In this Paper we obtain the analytical solution for the wave equation with the periodical dielectric 
function. We also calculate the frequency corresponding to the electromagnetic wave localized due to a defect in the 
array of the stacks of graphene stripes separated by a dielectric stripes. 

The Paper is organized in the following way. In Sec. [IT] we obtain the photonic band structure of the ID graphene- 
based photonic crystal. In Sec. Mil we find the frequency corresponding to the electromagnetic wave localized by the 
defect of a ID graphene-based photonic crystal. Finally, the discussion of the results and conclusions follow in Sec. lIVI 
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FIG. 1: The ID photonic crystal with graphene stripes, (a) the top view; (b) the side view, 
stripes between graphene stripes can be the same as the material of the dielectric substrate. 
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II. THE WAVE EQUATION FOR ID PHOTONIC CRYSTAL WITH GRAPHENE STRIPES 



We consider polarized electromagnetic waves with the electric field E perpendicular to the plane of graphene stripes. 
The wave equation for the electric field in a dielectric media has the form 34] 



AE(r. t) 



r(r,i) <9 2 E(r,i) 
~ 2 dt 2 







(1) 



where e(r, t) is the dielectric constant of the media, and c is the speed of light in vacuum. Looking for solutions with 
harmonic time variation of the electric field, i.e., E(r, t) — E(r)e lwt , and considering the propagation of wave in the 
a;— direction along the plane of graphene stripes and perpendicular to the graphene-dielectric boundaries one obtains 
from Eq. (P 



d 2 E z (x) uj 2 s(x,uj) 



dx 2 c 2 

The dielectric constant of the ID periodic structure is given by 



E z (x) = 



(2) 



e(x, uj) 



£o, for — |(a — b) + na < x < |(a — b) + na , 
Si(ui), for | (a — b) + na < x < i(a + b) + na, 



where £o is the dielectric constant of the dielectric, £i(lu) is the dielectric function of graphene multilayers separated 
by the dielectric barriers, a is the period of ID array of graphene stripes, b is the width of graphene stripes, and n is 
an integer. The ID photonic crystal with graphene stripes separated by dielectric layers with a thickness d is shown 
in Fig. [TJ By introducing the filling factor / the relation between a and b can be written as b = af. 

The dielectric function Sx(u) of a graphene multilayers system separated by dielectric layers with a dielectric 
constant £q and a thickness d is given by |35l l36j 



£i(w) = e + 



ud 



(3) 



where u g (uj) is the dynamical conductivity of doped graphene for the high frequencies (oj 3> kvF, w ^> t 1 ) at 
temperature T given by [H, HB| 
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Here e is the charge of an electron, r _1 is the electron collision rate, ij is the Heaviside step function, ks is the 
Boltzmann constant, and /i is the the chemical potential determined by the electron concentration, which is controlled 
by the doping. The chemical potential can be calculated as fx = (nno^^hvF, where the electron concentration is 
given by no — (/j,/(hvp)) 2 /n and vf — 10 8 cm/s is the Fermi velocity of electrons in graphene [36j . 

We solve the wave equation ([2]) to find the eigenfrequencies corresponding to the electromagnetic wave penetrating 
in the photonic crystal shown in Fig. [T] This wave equation is mathematically similar to the Schrodinger equation for 
an electron moving in a one-dimensional rectangular periodic potential barrier described by the Kronig- Penney model. 
The eigenenergies of the Schrodinger equation corresponding to the Kronig-Penney model presented in Eq. (|A1[) are 
given by [37j : 



cos(fca) = cosh(/36) cos [aa(l — /)] + 



a 



/3 2 



2af3 



sinh(/3&) sin [aa(l — /)] 
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where the wave vector k is in the range < k < 2n/a, and a and /3 are defined as 



UJ 



Applying Eq. ©, Eq. can be written in a form 
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The eigenfrequencies of the ID photonic crystal as functions of the wave vector k can be obtained by substituting the 
dielectric constant of the multilayer graphene given by Eq. (3) into Eq. ([7]). The solutions of Eq. ((TJ) for uj as functions 
of k provides the frequency band structure for ID photonic crystal formed by the periodically located parallel stacks of 
alternating graphene and dielectric stripes embedded into a background dielectric medium and allow one to determine 
the photonic band gap f2. Eq. ([7]) shows that the photonic band gap depends on the filling factor / and the period 
of ID array of graphene stripes a, as well as on the thickness d of the dielectric stripes that separate the graphene 
stripes. Let us mention that dependence of the photonic band gap on the thickness d of the dielectric layer, which 
separates the graphene stripes, is presented in ([7]), since the dielectric function £\{uj) ([3]) of a graphene multilayers 
system separated by dielectric layers depends on d. 

III. LOCALIZATION OF THE ELECTROMAGNETIC WAVE ON THE DEFECT 



Let us consider a defect in the array of stacks of alternating graphene and dielectric stripes embedded into a 
background dielectric medium. This defect is formed by one empty space or "ID vacancy" due to the absence of a 
stack of alternating graphene and dielectric stripes in one place, where it should be placed due to the periodicity. 
This place is filled by the dielectric. This extra dielectric stipe contributes to the dielectric contrast that results by 
adding the term — uj 2 /c 2 (si(uj) — eo)j(\x — xo\)E z (x) to the r.h.s. in Eq. @. Here 7(|a; — Xo\) = 1 for \x — x \ < b, 
and j(\x — xq\) — for \x — xo\ > b, where Xo corresponds to the coordinate in the middle of ID defect, which is the 
coordinate of the middle of the absent graphene stripe. As a result we obtain the wave equation for the electric field 
for ID photonic crystal with the defect: 



d 2 E z (x) uj 



dx 2 



+ — (e(x,uj) - (ei(w) - £obO - x \)) E z (x) = 



(8) 



Eq. ([5]) describes the periodic array of stacks of alternating graphene and dielectric stripes with the defect formed by 
one stack of alternating graphene and dielectric stripes being absent. 

We solve the wave equation ((8]) to find the eigenfrequency corresponding to the electromagnetic wave localized at 
the defect formed by a background dielectric medium due to the absence of the stack of alternating graphene and 
dielectric stripes in one place, where it should be placed due to the periodicity. This wave equation is similar to the 
Schrodinger equation describing the electron in the periodic potential energy in the presence of the potential energy 
of the defect placed at the point xq given in Appendix B. From the mapping of Klein-Gordon type equation given by 
Eq. (|B6[) the wave equation for the electric field in graphene-based photonic crystal has the following form 



4fc2 c 4 d 2 E z {x) 
"~3 dx 2 



2uj 2 n 2 ( £l (uj) - s h(\x - x Q \)E z (x) = (uj a - n 4 )E z (x) 
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where fco = 27r/a is the vector of the ID reciprocal lattice. This Klein-Gordon type equation has the eigenvalue 
w 4 — f2 4 . In Eq. f2 is the width of the forbidden band (photonic gap) in the spectrum of the electromagnetic wave. 

The electric field can be obtained by mapping Eqs. (|B7|) and (|B9|) that correspond to wavefunctions of the stationary 
states. The continuity of the electric field and its derivative in the points x — xq + b/2 and x = xq — b/2 result in the 
transcendental equation determining the spectrum of the even states: 

v/w 2 n 2 (ei(cd) - £ ) 2 - w 4 + ft 4 tan (V^ 2 (£iM - £q) 2 - ^ + ft 4 & 2 /c 2 ) = v 7 ^ 4 - ^ 4 | , (10) 

as well as of the odd states: 

v /a; 2 2 (e 1 (a;) - e ) 2 - w 4 + ft 4 cot - e ) 2 - w 4 + fi 4 6 2 /c 2 ) = -Vl^ 4 -^ 4 | • (11) 

Solving Eqs. (|10[) and (|lip with respect to the frequency w, we obtain the frequency corresponding to the photonic 
mode localized by the defect. 

IV. DISCUSSION 

In the calculations below we assume uq — 10 11 cm~ 2 . For simplicity, we consider the same material for a background 
dielectric medium and dielectric stripes between the graphene stripes. As the dielectric material we consider SiC>2 
with the dielectric constant £o = 4.5. Using Eq. ([7]) we calculated the band structure for ID graphene-based photonic 
crystal. The results of the calculations of the dispersion relation of the photonic crystal are presented in Fig. [2] 
The photonic band structure is calculated for different distances between graphene layers d. In our calculations we 
used the chemical potential for the electrons in graphene jU = 3.525 x 10~ 21 J, temperature T = 300 K, the period 
of the ID graphene stripes array a = 25 x 10~ 6 m, and the filling factor / = 0.3927. According to the results of 
our calculations, the photonic band structure almost does not depend on Eq due to the fact that £o <C |ei|. The 
results of our calculations demonstrate the strong dependence of the photonic band structure on the thickness d of 
the dielectric that separates the graphene stripes. At d = 1 nm and d = 5 nm the distance between the lower and the 
upper dispersion curves is larger in the middle than at the edges. At d = 3 nm and d = 10 nm the distance between 
the lower and the upper dispersion curves is larger at the edges than in the middle. 

Using Eqs. (flU]) and (ITU we calculate frequencies corresponding to the electromagnetic wave localized by the defect 
in the photonic crystal formed due to the absence of a stack of alternating graphene and dielectric stripes in one 
place, where it should be placed due to the periodicity. For the photonic crystal with the defect at d = 1 nm we 
have v = 1.79 THz, at d = 3 nm we have v = 1.85 THz, at d — 5 nm we have v = 3.91 THz, at d = 10 nm we 
have v = 3.62 THz. All these frequencies are located inside the photonic band gap. Let us mention that since our 
approach is based on the Luttinger-Kohn model [39l [ioj , the localized frequency v does not depend on momentum. 
The frequency localized by the defect is just a constant in the photonic band structure, which can be controlled by 
the thickness d of the dielectric stripes that separate the graphene stripes, as well as by the filling factor / and the 
period of ID array of graphene stripes a. 

In conclusion, the graphene-based photonic crystal proposed in this paper is a novel type one-dimensional photonic 
crystal formed by the array of the stacks of alternating graphene and dielectric stripes embedded into a background 
dielectric medium. This system can be analyzed as a ID photonic crystal and the corresponding wave equation for the 
electromagnetic wave propagating in such a structure can be solved in the framework of the Kronig-Penney model. 
The frequency band structure is obtained analytically as a function of the filling factor, the period of ID array of 
graphene stripes and the thickness of the dielectric between the graphene stripes. We obtain the photonic frequency 
corresponding to the electromagnetic wave localized at the defect of a photonic crystal formed by a background 
dielectric medium due to the absence of the stack of alternating graphene and dielectric stripes in one place, where it 
should be placed due to the periodicity. 
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Appendix A: The Schrodinger equation with the periodic potential 

The wave equation @ with the dielectric constant given by Eq. Q can be mapped onto the Schrodinger equation 
for an electron in the ID periodic potential function approximated by a rectangular potential barrier, used in the 
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FIG. 2: The dispersion relation for the ID graphene-based photonic crystal with a — 25 fim and filling factor / = 0.3927 for 
different thicknesses of the dielectric between graphene stripes. 



Kronig-Penney model [37[ : 



h 2 d 2 ^{x) _,, , . 1, 1, 

- Eip(x) , for (a - 6) + na < x < -(a - 6) 



2m <9x 2 w ' 2 V ; 2 

f 2 



na. 



d 2 tb(x) ?i 2 u; 2 1 1 

1 V(x) = Etplx) , V(x) = for -(a - 6) + na < x < -(a + b) + na, (Al) 



2mo 9x 2 ' 2moc 2 

where mo is the mass of an electron, ip(x) and E are the wave function and energy of an electron, V(x) is the potential 
due to the ion of the ID crystal lattice approximated by a rectangular barrier, a is the period of ID array of the 
scatterers with the rectangular potential, b is the width of the rectangular potential barrier, and n is an integer. 

It is easy to see that the wave equation with the dielectric constant given by Eq. ([3]) can be mapped onto the 
Schrodinger equation (|A1|) for an electron in the ID periodic potential corresponding to the Kronig-Penney model. 
The mapping relations are the following 



r/>(x) = E z (x) , 

a = —2~E = —uj , A2) 



—2-(E-V{x)) EE — —I 



The solution of Eq. (IA1[) is well known and is given by expression ([5]) . 
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Appendix B: The Dirac-type equations for the the electromagnetic wave in the photonic crystal with the 

defect 



Following Ref. [l5j we can map the wave equation for the photonic crystal with the defect onto the Schrodinger 
equation in the periodic field with the defect. After mapping of Eq. (|5J) onto the Schrodinger equation describing the 
electron with the effective electron mass mo in the periodic potential energy V(x) in the presence of the potential 
energy of the defect W7(|x — xq\) placed at the point xq we have 



h 2 d 2 

2mo dx 2 



In Eq.dEHl ip(x)=E z (x), 



V{x)-Wj{\x-x \) 



h 2 u 2 



(Bl) 



and the potential V(x) is given in Eq. (Al) and potential W is defined as 



W = 



2m c 2 



(ei(w) - £ ) 



(B2) 



Eq. (|B1[) has the same form as Eq. (6) in Ref. (38[. However, in our case the potential W is defined by Eq. 
and corresponds to the potential of the defect in the Schrodinger equation for an "electron" in the periodic field of 
the crystal lattice and in the presence of the defect. 

We will reduce the problem of the Schrodinger equation for a particle in the periodic potential V(x) related to the 
system of the periodically placed stacks of alternating graphene and dielectric stripes embedded into a background 
dielectric medium and an "defect potential" W related to a defect to a much simpler equation for the envelope 
wavefunctions. Taking into account the two-band structure and following the two-band model f38j ] we introduce a 
two-component spinor wave function ijj(x) 



i/j(x) = 



tp(x) 



(B3) 



where two different neighboring bands are described by wave functions <p(x) and x( x )- Note that Eq. (|Blj) contains 
both the periodic function V(x) corresponding to the ideal lattice and W, describing the potential of a defect. Without 
a defect the energy spectrum would be described by two neighboring bands and the gap between them. Applying the 
standard two-band approach, we obtain an effective Dirac-type equation for the envelope spinor wave function, which 
implies the periodicity provided by V(x) [l5ll38j. 

Note that Eq. (|B1[) describes an electron in the periodic potential of the ideal crystal lattice V(x) and the potential 
of the defect W. If the solution corresponding to the absence of impurity W — is known, the energy levels of 
the electron localized by the defect can be obtained by replacing Eq. (jBll) by the Dirac-type equation according to 
Luttinger-Kohn model described in Rcfs. [38l - l40j . This model implies a Dirac-type equation for the two-component 
spinor wave function. According to Ref. |38|. the functions (f n (x) and Xn(x) defined as 



<Pn{x) = *^2c n (k)exp[ikx/h] , 
Xn{x) = dn{k') exp [ik'x/h] 



(B4) 



satisfy the set of the second order partial differential equations. Considering only two neighboring bands corresponding 
to the wave function tp(x) given by Eq. (|B3[) this set of equations for (p n (x) an d Y n (x) can be reduced to the Dirac-type 
equations for the two-component spinor (|B3I) . which has the following form [38J 



A w - Wj(\x - xq\)] tp(x) + ihsa- 
A w - Wj(\x - scq |)] x{x) + ifiso. 



dx(x) 

dx 
dip(x) 

dx 



, 
, 



(B5) 
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where s = hko/(V3mo). In Eq. (IB5|) a x is Pauli matrix and A w is the width of the forbidden band in the electron 
spectrum that, as it follows from the mapping of the wave equation for the electric field graphene-based photonic 
crystal onto Eq. (|B1|) . is A u — g^^a , where fl is the width of the forbidden band (photonic gap) in the spectrum 
of the electromagnetic wave. Eqs. (|B5|) are obtained in the limit — A 2 J \/(2A 2 J ) <C 1. Defining the effective 
mass of a quasiparticle as = 377igA [i ,/(?i 2 fcQ) , and following the standard procedure of quantum electrodynamics 
[4ll |42| we obtain from the system of Dirac-type Eqs. (|B5[) the following Klein-Gordon type equation. Note that 



the Klein Gordon equation can be reduced to a Schrddinger-like equation with an effective energy and an effective 
potential. If the potential is weak enough to ignore the W 2 term, the relativistic formalism becomes equivalent to the 
non-relativistic formalism. More importantly, in situations where the Klein-Gordon equation is not exactly solvable, 
the Schrodinger form of the Klein-Gordon equation sheds some light on the problem as it could be reduced to a 
solvable eigenvalue problem: 



-?iV— - 2m LJ s 2 W 1 (\x - x \) 
ax* 



(B6) 



This Klein-Gordon type equation has the form of the 2D Schrodinger equation for a particle in the square well 
potential well with the eigenvalue 



£, 



A 



2m^s 2 

The wavefunctions of the even stationary states have the form (£ u < 0): 



^(x) = Ai cos 



^(x) = Biexp 



yj2m u (W-£,„)/h 2 x 



\x — xq\ < b 



x \ > b 



(B7) 



The continuity of the wave function and its derivative in the points x = xq + b/2 and x = xq — b/2 result in the 
transcendental equation determining the spectrum of the even quantum states: 



y/W - £ u tan sJ2m u {W - £^)b 2 /fi 2 = y/\£i 
The wavefunctions of the odd stationary states have the form (£ u < 0): 



(B8) 



*(x) = A 2 sin 
*(ar) = i3 2 exp 



^2m^W -£„)/ <H 2 x 



—>J2mu£ u /Ti x 



\x — xq\ < b 
\x — xq\ > b . 



(B9) 



The continuity of the wave function and its derivative in the points x = xq + b/2 and x = xq — b/2 result in the 
transcendental equation determining the spectrum of the odd quantum states: 



VW-£ U cot fimJW -^W/h 2 
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